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NOTE SUR LA DETERMINATION DES PARAMETRES DU 
MODELE DE BURSA- WOLF 

ABDELMAJID BEN HADJ SALEM 



1. Introduction 

Avec le developpement de la technologie de positionnement spatial (GPS, 
GLONASS, Galileo, ComPass), laquelle fournit a I'utilisateur sa position 
(X, Y, Z) tridimensionnelle dans un systeme geocentrique mondial donne, 
par exemple pour la technologie GPS c'est le systeme dit WGS84: (World 
Geodetic System 1984), il est necessaire de savoir la transformation de pas- 
sage du systeme geodesique mondial an systeme geodesique national ou local. 
Nous presentons ci-apres en detail comment determiner manuellement les pa- 
rametres du modele de Bursa- Wolf de transformations de passage entre les 
systemes geodesiques. 

Nous utilisons par la suite les notations suivantes : 

- (Xi , Yi , Zi ) les coordonnees cartesiennes 3D dans le systeme geocen- 
trique (O, Xi, Yi, Zi) (systeme 1), 

- (X2,Y2, Z2) les coordonnees cartesiennes 3D dans le systeme local (sys- 
teme 2) {0',Xi,Yi,Zi). 

2. Le Modele de BURSA - WOLF 
Ce modele s'ecrit sous la forme vectorielle : 

(1) X2 = T + {l + m).R{rx,ry,rz).Xi 
ou : 

- X2 est le vecteur de composantes {X2,Y2, Z2Y' , T designe transposee, 

- T = O'O est le vecteur translation de composantes {Tx ,Ty ,Tz)'^ entre 
les systemes 1 et 2, 

- 1 + m est le facteur d'echelle entre les 2 systemes, 

- R{rx, ry, rz) est la matrice de rotation (3 x 3) pour passer du systeme 1 
au systeme 2, 

- Xi est le vecteur de composantes {Xi,Yi, Zi)'^ . 

En developpant ([T|), on obtient : 

/X2\ /Tx\ ( 1 -rx ry\ ( XA 

(2) Ya = Ty + (1 + m) rx 1 -rz Yi 
\Z2j \Tzj \-ry rz 1 / \Zij 

avec (rx, ry, rz) les rotations comptees positivement dans le sens contraire 
des aiguilles d'une montre. Comment determiner les parametres modele ([T]) ? 
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3. Determination de l'echelle 1 + m 

On suppose donne un ensemble de points Pi pour i = l,n connus dans 
les deux systemes et 82- On ecrit I'equation ([T]) pour deux points Pj et 
Pk, d'ou : 

(3) X{Pj)2 = T + {l+m).R{rx,ry,rz).X{Pj)i 

(4) X{Pk)2 = T + (1 + m).R{rx, ry, rz).X{Pk)i 
Par difference, on obtient : 

(5) (PjPk)2 = (1 + m).R{rx, ry, rz).(PjPk)i 

On prend la norme des deux membres de ^ et que 1 + m > : 
(6) 

||(PjPk)2|| = ||(l+m).i?(rx, ry,rz). (PjPk)ill = (l+m)||ii(rj;,ry,rz).(PjPk)i| 

Comme R est une matrice de rotation, done son application a un vecteur est 
une isometrie, c'est-a-dire qu'elle laisse invariant la norme ou la longueur du 
vecteur soit : 

(7) = ||X|| VX 
On a done : 

(8) ||(PjPk)2|| = (l+m)||(PjPk)i|| 
Soit : 

(9) l + .n = lfM^™ 

iV^||(PjPk)i|| 

N designe le nombre de couples de points PjPk,j ^ k. 

4. Determination des rotations {rx,ry,rz) 

Connaissant (1 + m), pour un couple de points Pj, Pk, on a : 

(PjPk)2 = {I + m).R{rx,ry,rz).{PjPk)i 

Detaillons la matrice R : 
(10) 

R = \ rx 1 -rz I = I 1 I + I rx -rz \ = I3 + Q 
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avec I3 la matrice Unite et Q la matrice : 








Alors I'equation ([5]) devient : 

(12) (PjPk)2 = (1 + m).(/3 + Q{rx, ry, rz)).(PjPk)i 
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Soit comme m « 1 et m << 1 : 

(13) Qirx, ry, rz).(PjPk)i = (1 - m).(PjPk)2 - (PjPk)i 

En posant : 
(14) 

/AXjA /AX, A /v^ = il-m)AX'^,-AX,k^ 

(PjPk)2 = Aigj ; (PjPk)i = Ay, J ; v=\v2 = {l- m)AY;, - AY,k 

\AZ'^J \AZ,kJ V^3 = (1 - m)AZ'^, - AZ,k, 

Alors, on obtient I'equation : 

(15) Q(rx, ry, rz).(PjPk)i = v 

ou encore : 




(16) 

\ -AX,fc 
Or le determinant de la matrice Q' : 

(17) Q' 



rx\ 




(vi 


ry 


\-\ 




rzj 




Us 




est nul. Pour passer de cette consequence, on utilise pour chaque ligne du 
systeme (I16p un couple de points ij ce qui donne le systeme : 



(18) 



Le systeme (llSp devient resolvable ce qui permet de determiner les trois 
rotations rx, ry et rz. 




rx\ 




( Vjkl 


ry 




Vlm2 


rz J 




Wm3 



5. Determination des composantes de la Translation T 

Les composantes Tx, Ty, Tz du vecteur translation sont determinees a 
partir des coordonnees des points Pj connus dans les deux systemes a partir 
de : 

(19) Txj = X2j - (1 + m){Xij - rxYij + ryZij) 

(20) Tyj =Y2j-{l+ m){rxXij + Yij - rzZ^j) 

(21) Tz^ = Z2j - (1 + m){-ryXij + rzYij + Zy) 
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Les composantes Tx, Ty, Tz sont obtenues par une moyenne sur les N points 
communs a savoir : 



(22) 


Tt — 


(23) 




(24) 
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